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Abstract 



Behaviour of a weekly self-interacting scalar field with a small mass in the de 
Sitter background is investigated using the stochastic approach (including the 
case of a double- well interaction potential). Existence of the de Sitter invari- 
ant equilibrium quantum state of the scalar field in the presence of interaction 
is shown for any sign of the mass term. The stochastic approach is further 
developed to produce a method of calculation of an arbitrary anomalously 
large correlation function of the scalar field in the de Sitter background, and 
expressions for the two-point correlation function in the equilibrium state, 
correlation time and spatial physical correlation radius are presented. The 
latter does not depend on time that implies that the characteristic size of do- 
mains with positive and negative values of the scalar field remains the same 
on average in the equilibrium state in spite of the expansion of the t = const 
hypersurface of the de Sitter space- time. 
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I. INTRODUCTION 

De Sitter space-time is one of the most fundamental and symmetric space-times. It is 
a constant-curvature space-time completely characterized by only one constant H } and its 
group of symmetries is as large as the Poincare group of symmetries of the Minkowski space- 
time (10 generators). That is why it is so significant to get a clear understanding of behaviour 
of quantum fields including interacting ones in this background. This program was launched 
in the pioneer paper by Chernikov and Tagirov [1] for the case of a free massive scalar held 
(see also [2]- [7] and other papers on this subject). Also, quantum held theory in the de 
Sitter background has extremely important cosmological applications in connection with 
the inflationary scenario of the early Universe and generation of scalar perturbations (both 
adiabatic and isocurvature ones) and gravitational waves from vacuum quantum fluctuations 
during an inflationary stage. 

Due to the high symmetry of the de Sitter background it is natural to define the back- 
ground (or "vacuum") state of any quantum field as a state vector invariant under action 
of the full group of symmetries of the de Sitter space-time 0(4, 1). Note that, due to the 
possibility of a static representation of a part of the de Sitter space-time surrounded by 
the event horizon, the same state, if exists, may be also called a "thermal" state with the 
temperature T = H/2tt [8]. This fact, however, simply reflects the symmetries of the de 
Sitter space-time and does not mean that average values of operators in this state, e.g. (</> 2 ), 
are given by corresponding thermal values in the flat space-time with the Gibbons- Hawking 
temperature (in fact, the former are much larger in the case we are interested in). One may 
also expect that this background state is a stable attractor for an open set of other quantum 
states in the Schrodinger picture, in other words, that it is a stable equilibrium state. In 
the Heisenberg picture, this means that average values of all possible operators in an open 
set of constant state vectors approach those in the equilibrium state as \t\ — > oo. 

The investigation of a free massive scalar field in the de Sitter background has shown 
that: 

1) such an equilibrium background state does exist if m 2 > 0, where m 2 is the effective mass 
squared in the de Sitter background including a contribution from the non-minimal coupling 
term £R<f> in the dynamical equation for the scalar field (if such term is present), R being 
the scalar Ricci curvature; 

2) there is no de Sitter invariant state for m 2 < 0; 

3) in the case ra = 0, if the de Sitter (inflationary) stage begins (or is "switched on") at the 
moment t = , the average value (<f> 2 ) in the Heisenberg state which was the vacuum state 
for t < (according to, say, the quasi-adiabatic definition of vacuum) and in an open set of 
close states grows as H 3 t/4:Tr 2 for Ht >> 1 [4-6]; 

4) for m 2 > 0, m <C H, the value of (<f> 2 ) in the de Sitter-invariant state is anomalously large 
as compared to the corresponding thermal value with the Gibbons-Hawking temperature, 
namely, (<p 2 ) = 3H 4 /8ir 2 m 2 [3]. 

Now the question arises what happens in the case of a self-interacting scalar field in the 
de Sitter background. The most natural type of such an interaction is the quartic coupling 
A</> 4 /4 with < A <C 1 in the Lagrangian density of the scalar field. This problem was 
first considered in [2]. More detailed recent investigations [9,10] showed that there exists 
a breakdown in the perturbation expansion around the free-field (A = 0) operators due to 



infrared divergences if m 2 < 27H 2 /1Q. This shows that the held theory has a completely 
non-trivial infrared structure for these values of m 2 (especially in the most interesting case 
\m 2 \ <C H 2 ) that lies beyond any finite order of the perturbation expansion in A. 

Fortunately, there exists a method to manage this difficulty and to solve the problem in 
the case \m 2 \ <C H 2 that uses specific properties of the de Sitter space-time (so it cannot 
help infrared problems in the flat space-time). This is the stochastic approach to inflation 
(also called "stochastic inflation") that is based on the idea first proposed in [5] that the 
infrared part of the scalar field may be considered as a classical (i.e, c-number) space- 
dependent stochastic field which satisfies a local (i.e., not containing spatial derivatives) 
Langevin-like equation. The final form of the corresponding Fokker-Planck equation for 
the probability distribution function (PDF) of this classical field which is suitable for our 
purposes was derived in [11,12] (a partial form of the latter equation for the case of a 
massless non-interacting scalar field was presented in [13]). In this paper, we make a further 
step and expand this method to the calculation of two-point and higher-point correlation 
functions that enables to find any anomalously large Green function of the scalar field in the 
equilibrium background quantum state (large as compared to the corresponding thermal 
quantity in the flat space-time). Such functions are generally average values of arbitrary 
functions of the scalar field itself taken at different points in the de Sitter space-time. On 
the other hand, any non-zero average values including second time derivatives, squares of 
first time derivatives and spatial gradients of the field are not anomalously large and cannot 
be exactly calculated by the method used. 

So, the structure of the paper as follows. In Sec. 2, we present basic field equations and 
make a calculation of the equilibrium value of (<f> 2 ) using the Hartree-Fock (or Gaussian) ap- 
proximation to get a crude estimation of this quantity. In Sec. 3, we introduce the stochastic 
approach, display the Fokker-Planck equation and find a static equilibrium solution for the 
one-point PDF. Expressions for the temporal autocorrelation function of the scalar field in 
different cases are derived in Sec. 4. In Sec. 5, the general two-point function is found and 
its de Sitter invariance in the equilibrium state is proved. Sec. 6 contains discussion of a 
general picture of fluctuations of the self-interacting scalar field in the de Sitter background. 

II. BASIC EQUATIONS AND THE SCALAR FIELD EVOLUTION IN THE 
HARTREE-FOCK APPROXIMATION 

The Lagrangian density of a massive scalar field with the quartic self-interaction is (c = 
fi = 1 and the Landau-Lifshits sign conventions are used throughout): 

C = \(^^-m 2 J 2 + m 2 )-^-. (1) 

In the de Sitter background, the scalar curvature R = — 12H 2 } H = const, so the effective 
mass in the limit A = is m 2 = m 2 , + 12^H 2 . Thus, the wave equation for the scalar field 
in this background is 

<^ + m 2 ^ + A<^ 3 = (2) 

which is also the equation for the corresponding Heisenberg field operator, here ; denotes the 
covariant derivative. We assume further that \m 2 \ <C H 2 . Properties of the quantum field 



theory with the Lagrangian density (1) at energies E 3> H essentially coincide with those 
in the flat space-time because corresponding modes of the held do not feel the space-time 
curvature. So, this theory is renormalizable but it possesses the Landau "zero-charge" pole 
at very high energies. Thus, it may exist (both in the flat and curved space-time) as a part 
of a more fundamental underlying theory only. Also, A is running with energy E for E 3> H 
according to the flat space-time renormalization-group equation. The de Sitter curvature H 
produces an effective infrared cut-off for all logarithmic radiation corrections (see, e.g., [5]), 
so by A we shall further mean X(H). To avoid problems with large radiation corrections we 
assume that < X(H) <C 1 (then A(0) is even less) and that the energy H, though possibly 
being high and even not far from the Planckian one Epi = Mpi = 1/y/G, is much lower than 
the energy of the Landau pole that is exponentially large. 
We take the metric of the de Sitter space-time in the form 

ds 2 = dt 2 — a^e 2 f o?x 2 = (o??y 2 — o?x 2 ) , a = const (3) 

[Hi]) 2 

where rj is the conformal time defined by a e Ht = l/(—Hr]) for — oo < rj < 0. This metric 
covers only one half of the whole de Sitter space-time. This is, however, sufficient for our 
purpose because if a de Sitter invariant state is constructed in a part of the space-time, it 
may be analytically continued to the whole space-time. 

As a first step to the understanding of the behaviour of the interacting scalar field (1) in 
the de Sitter background, let us calculate the evolution of (<f> 2 ) in the Hartree-Fock approxi- 
mation following [4,5] (for a recent discussion of this method, see [14]). This approximation 
is also called Gaussian because the wave functional in the Schrodinger representation is as- 
sumed to be Gaussian. For us, however, it will be more convenient to work in the Heisenberg 
representation. The Hartree-Fock approximation is known to be exact in the non-interacting 
case (A = 0). Also, it becomes exact if we consider a multiplet of N scalar fields with the 
O(N) symmetry and take the limit N — > oo (the so-called 1/A expansion). For A ^ 
and N = 1, one may expect to get a qualitatively correct result only (if at all). Still the 
Hartree-Fock approximation is so widely used that it is instructive to see what it produces 
in our case. 

To find the evolution equation for (</> 2 ), let us multiply Eq. (2) (considered as the operator 
equation) by <f> and then take its average in some state vector. We get 

\{m = ^^)-m 2 ^ 2 )-X{^). (4) 

It is known from the analysis of the free massless case [4-6] that the main contribution to 
((f) 2 ) is made by long-wave, slowly changing with time modes. So, only the term 3i7^ may 
be kept from the d'Alambertian operator in the left side of Eq. (4). The first term in 
the right side of Eq. (4) contains two derivatives and does not have an anomalously large 
infrared part. Let us use the Hartree-Fock (or Gaussian) approximation to estimate the last 
term in Eq. (4). If we assume that (</>) = initially, it will be formally zero at all times due 
to the reflection symmetry <f> — > —<f>. Then we get (</> 4 ) = 3(</> 2 ) 2 as if the operator <f> would 
be a Gaussian stochastic quantity. So, Eq. (4) reduces to 

f) H 3 2m 2 2X 



where the first term in the right side is introduced to recover the correct linear growth of 
(<f> 2 ) in the limit A — > 0. It originates from the time dependence of the upper cut-off for 
the large infrared part of (<f> 2 ) in the space of conformal wave vectors k : k ~ aH , that 
corresponds to a time-independent cut-off in the physical space. On the other hand, the 
lower cut-off is constant in the &-space because it is determined by the moment when the 
de Sitter expansion begins (cf. [4,5]). It is clear that all solutions of Eq. (5) approach a 
constant equilibrium value at t — > oo. In particular, ((f) 2 ) —■ H 2 /irV8\ ~ 0.113 H 2 / y A if 



m 2 = 0. This corresponds to the appearance of the effective mass rn 2 jf ~ H 2 y\ due to 
self-interaction. 

Certainly, we cannot expect that the Hartree-Fock approximation, being exact for free 
quantum fields, produces a quantitatively correct result in the presence of self-interaction. 
In particular, the form of the PDF of <f> that it assumes is completely misleading, it is not 
Gaussian at all. So, let us turn to the quantitatively correct treatment of the problem now. 



III. STOCHASTIC APPROACH AND THE FOKKER-PLANCK EQUATION 

In the stochastic approach [11,12] we represent the Heisenberg operator of the quantum 
held 6 as 



x, t) = <Kx, t) + / -^j6(k - ea(t)H) W^(t)e~^ + 4^(<)e«' kx l . (6) 



Here </>(x, t) is a coarse-grained or a long wavelength part of </>, 6(z) is the step function 
and e is a small constant parameter which in fact is not arbitrarily small (we shall return 
to this point below and determine the inequality it should satisfy). Note that if we consider 
a quasi-de Sitter background with a slowly changing curvature \H\ <C H 2 } then eH should 
be taken as a constant to avoid the appearance of additional (though small) terms in the 
Langevin and the Fokker-Planck equations written below. So, <f> is the held <f> averaged over 
a constant physical 3D volume slightly larger than the volume inside the event horizon. 

The short wavelength counterpart is characterized by the mode functions <f>^(t) with a^ 
and 4 being annihilation and creation operators, respectively. It is small, so it satisfies the 
linear massless equation 

4 + 3#4 + -—-M = 0, (7) 



a 2 (t) 



which is solved as 



\H{-tfHf{-kn) = ^(v~l) e~ tk \ (8) 



?k = 
with Hi (z) being the Hankel function of the hrst class. With the above choice of mode 

2 

functions, the constant Heisenberg quantum state annihilated by a^'s corresponds to the 
usual adiabatic (or Minkowski) vacuum in the limit rj — > — oo. 

From the equation of motion of the scalar held with some interaction potential V(<f>), 
□ </> + V'(cf>) = 0, we obtain the following equation for the slowly- varying coarse-grained part 



i(x,*) = -g^n^)+/(x,<), (9) 

where /(x, t) is given by 

/(x,*) = ea(t)H 2 f -^j6(k - ea(t)H) W^(t)e~^ + a{f k (t)e tkx ] . (10) 

J I/7TI2 L J 



Eq. (9) may be regarded as the Langevin equation for the stochastic quantity <f> (the fact 
that it is still an operator becomes irrelevant) with a stochastic noise term /(x, t) whose 
correlation properties are given by 

±± sin z 
(/(x 1 ,t 1 )/(x 2 ,t 2 )) = — S(h -t 2 )j {ta{t)H\^ -x 2 |), j (z) = . (11) 

471"^ z 

We note that the white-noise property of the above correlation is preserved even if we take 
into account the first-order correction arising from self-interaction. 

It is straightforward to derive the Fokker-Planck equation for the one-point (or one- 
domain) PDF pi[<f)(x.,t) = ip] = pi[ip(x.,t)] = pi((p } t). It has the form: 

d/?i[y?(x,*)] 1 d ff 3 dVi[y(x,t)] _ 

dt = 3~Hd^ {V MX'WpMxM + ^^ = r,PM^)]. (12) 

Once we have found a solution of the above equation, we can calculate the expectation value 
of any function of </>(x, t), (F[<f)(x.,t)]), using the expression 

(Ftf(x,t)]) = j d<pF(<p) Pl [<p(x,t)]- (13) 

Eq. (13) explains the sense in which the transition from quantum to classical behaviour 
takes place here: though the scalar held <f> remains a quantum operator formally, we can 
introduce an auxiliary classical stochastic scalar held <p(x.,t) with the PDF determined from 
Eq. (12) that has the same expectation values for all observables with excellent accuracy 
(the error is due to the decaying mode of <f> that is exponentially small, oc exp(—3Ht)). 

If H is not exactly constant, but slowly varying, it is more natural to write the Fokker- 
Planck equation in terms of the independent variable In a = / Hit) dt (see the discussion in 
[12]). Then the PDF P p ((p } lna) introduced as in [15] which refers to the "whole Universe" 
reduces to pi(ip,\na) after being normalized to unity according to the prescription given 
in [16]. Also, the distribution of domains N(ip, hi a) used in [17] is equal to a 3 / 0i(( J c,ln a), if 
one takes eH = const as pointed above. Clearly, in our case there is no difference between 
Ht and In a at all. 

In general, the solution of (12) can be written in the form 

Pifat) = exp (-^J^) E^Me-^', (14) 

where <& n (<p) is the complete orthonormal set of eigenfunctions of the Schrodinger-type 
equation 



1 d 



,2 



2 dip 2 
with the effective potential 



+ %) 



*»M = ^rMv). (15) 



^M = ^ [»'M 2 - »»] , v(<p) = ^V(<p). (16) 

The coefficients a n are given by an initial condition of p x at t = t as 

a n = Jd<p Pl (<p,tyM$ n (<p). (17) 

Because the left-hand-side of (15) can be recast in the form 

the eigenvalues A n 's are nonnegative. If 

/oo 
e- 2v ^d<p, (18) 

-oo 

is finite, we find A = and the corresponding eigenfunction is given by 

$ ( v? ) = JV-2 e -"H (19) 

Then the normalization condition f^dip pi(ip,t) = 1 results in a = N~? } so that (14) 
reads 

oo 

Pl {<p,i) = p lq (p) + e~ v M J2 «^.(^)e- A " (f - fo) - (20) 

71=1 

Therefore, any solution (14) asymptotically approaches the static equilibrium solution 

KM = iV" 1 exp ("|^m) = N-'e-W (21) 

This PDF is significantly non-Gaussian in the presence of self-interaction. In particular, if 
V(<p) = \p\ then 

, , /32vr 2 A\^ 1 / 2tt 2 A^ 4 \ , , 



One of distinctive non-Gaussian properties of the PDF (22) is its negative kurtosis K 
3 ~ —0.812. As for the dispersion of cp, it gives 



( V 2)2 



/ 3 T(1)H 2 H 2 , , 

Comparing it with the Hartree-Fock result for this quantity presented at the end of the 
previous section, we see that the latter is ss 15% less — not a bad result for the Hartree-Fock 
approximation. Certainly, it cannot be used to calculate such things as the kurtosis at all. 
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IV. EQUILIBRIUM TEMPORAL AUTOCORRELATION FUNCTION IN THE 

STOCHASTIC APPROACH 

Of course, the one-point PDF is not enough to obtain the whole picture of the scalar held 
distribution in the de Sitter space-time. We need to know all two-point and higher-point 
correlation functions. We shall show now that no new equations are necessary for this aim 
and that all these functions can be constructed from various solutions of the Fokker-Planck 
equation (12) (note that we are speaking about functions having anomalously large infrared 
part only). Let us begin with the temporal autocorrelation function of the scalar held. 

The autocorrelation function (or the temporal two-point function) G(ti 7 t 2 ) = 
(</>(x, ti)</>(x, t 2 )) can be calculated from the two-point PDF at equal space points 
p 2 [(f>(x,ti) = (fi,(f>(x,t 2 ) = <p 2 ] = /02[^i(x,ti),(^ 2 (x,t 2 )] as 



C(ti,t 2 ) = / dcp 1 dcp 2 cp 1 cp 2 p 2 [cp 1 (x,t 1 ),cp 2 (x,t 2 )]. (24) 

Here p 2 [<fi(x, t t ), ^ 2 (x, t 2 )] is given by 

/ o 2 [(^i(x,t 1 ),(^ 2 (x,t 2 )] = n[^ 1 (x,t 1 )|(^ 2 (x,t 2 )] / o 1 [(^ 2 (x,t 2 )]6'(t 1 -t 2 ) 

+Il[<p 2 (x,t 2 )\<p 1 (x,t 1 )]p 1 [<p 1 (x,t 1 )]9(t 2 - h), (25) 

where II[(^i(x, ti)\tp 2 (x. } t 2 )] is the conditional probability to find </>(x, ti) = ipi provided that 
</>(x, t 2 ) = (p 2 . It satishes the Fokker-Planck equation (12) with respect to both pairs of its 
arguments <p-\_,ti and (fi 2} t 2 [18]: 

dU, 

and 

dU, 



-[(^ 1 (x,t 1 )|^ 2 (x,t 2 )] = r vi n[^ 1 (x,t 1 )|^(x,t 2 )] (26) 



dt -[^i(x,ti)|^ 2 (x,t 2 )] = r V2 n[^ 1 (x,t 1 )|^ 2 (x,t 2 )], (27) 

with the initial condition 

n[(^ 1 (x,t 1 )|(^ 2 (x,t 1 )] = %! - Lp 2 ). (28) 

Since Y [fi is independent of time, we find il[( J c 1 (x, t 1 )|( J c 2 (x, t 2 )] depends on time only through 
l^i ~~ t 2 \ using (26) — (28). For t\ > t 2} / o 2 [( J c 1 (x,t 1 ),( J c 2 (x, t 2 )] itself satishes 

T^[^i(x,ti),(^ 2 (x,t 2 )] =T Vl p 2 [ip 1 (x,t 1 ),ip 2 (x,t 2 )]. (29) 

For a de Sitter invariant state, (</>(x, ti)</>(x, t 2 )) depends only on \t\ — t 2 |, so that 

G(h,t 2 ) = G(\h -t 2 \) = (?) C(\ti - h\), £(0) = 1 . (30) 

This is realized in (24) if and only if / o 1 [( J c(x, t)] is time- independent. In other words, in 
the language of stochastic inflation, the absence of a de Sitter invariant state for a free 
massless minimally- coupled scalar field is due to the non-existence of any static solution to 
(12) with V(cf>) = 0. On the contrary, one may expect to recover the de Sitter invariance of 
the two-point correlation function by adding the self-interaction V((f)) = ^</> 4 for which we 
have found the static equilibrium solution (22). 
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A. Massive non-interacting case 

Before proceeding to the explicit calculation of the autocorrelation function for the mass- 
less ^</> 4 theory, let us demonstrate that we can obtain the familiar formula for the Green 
function of a free massive scalar held by taking V((f)) = |m 2 </> 2 with < m 2 <C H 2 . Multi- 
plying (29) by (p\(p 2 an d integrating over ipi and tp 2} we find 

^mmh)) = -^(vwiMh)) = ~«ti)^)>- (3i) 

— 2 

Since the equilibrium value of {</> (£)) can be calculated as 

(32) 







<A*)> = 


= / dipip 2 


''PlM 


3H 4 

9>TT 2 m 2 


i 


. (31) is 


solved to 


yield 














(#i)?(i 2 )) = 


3H 4 

9>TT 2 m 2 


exp - 


m 2 


t 2 



(33) 

This has to be compared with the two-point Green function in the de Sitter invariant 
vacuum calculated by methods of the held theory [1-3], namely, 



H 2 (l-c)(2-c) ( l + z\ 3 /9 m 2 m 2 

(<£(xi,*iU(x 2 ,* 2 )) = t rF lc,3-c,2;^—) : c= - - \ = -, 

m i, i;n 2, 2,1 1Q7rs - m7r ^ _ c ) V' " 2 /' 2 V4 # 2 3# 2 ' 

(34) 

where F is the hyper-geometric function and z = z(x i} x 2 ) is the de Sitter invariant function 
of (xi,ti) and (x 2 ,t 2 ) dehned by 

H 2 

z(x 1 ,x 2 ) = coshi7(t 1 - t 2 ) age fftl+fft2 |xi - x 2 | 2 . (35) 

It is related to the geodesic interval s between (xx,^) and (x 2 ,t 2 ) by the formula z = 
1 + H 2 s 2 /2. For X! = x 2 and \t 1 — t 2 \ ~^> H _1 , the leading term of (34) reads 

Wl)< ^ 2)) ~ 8^ exp ("S^ 1 ~ hl ) ' (36) 

in agreement with (33). Since the autocorrelation function G(|t|) is the integral over ipi and 
(p 2 of some solution of the Fokker-Planck equation (12) (see (24)), it itself may be expanded 
using the eigenvalues of Eq. (15): 



G(\t\) = £ b n e~^\ J2 b n = (<A (37) 



-A„|t| 
n=l n=l 



(the coefficient b = due to the reflection symmetry). If A = 0, then A n = m 2 n/3H, n = 
0, 1, 2.... We see that in the absence of self-interaction the autocorrelation function contains 
a contribution from the only one eigenvalue Ai. Let us introduce the correlation time t c 
using the definition £(t c ) = 0.5. Then t c = In 2 A^ 1 = 3 In 2 Hm~ 2 > H' 1 for A = 0. 
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B. Massless interacting case 

Now we turn to the massless minimally-coupled A</> 4 /4 theory for which Eq. (31) is 
not readily soluble, and we must deal with the two-point PDF directly. Let us define the 
following dimensionless quantities 



Ht, (p=(^L^\ ^- = Q^ ^ = ^-^1, n^Q" 1 !!, h = Q~ 2 P2 (38) 



V 24tt 2 ' r V 3 J H 
with which the Fokker-Planck equation reads 

The corresponding dimensionless effective potential W((p) = Q~ 2 W(ip), see (16), is depicted 
in Fig. 1. 

Then the dimensionless autocorrelation function g(r) = (^(t + t)(/?(t )) = Q 2 G(t) is 
given by 

g(r) = / J^i^ 2 ^in[^l(To + T)|^ 2 (To)]/Oleq(^2)^2- (40) 

Using the Fokker-Planck equation for II and the initial condition n[(^ 1 (T )|(^2( T o)] = fi(<Pi ~~ 
(^2)5 we find 



d n g(r) 



dfadfcipiTl Il[^i(To)|^2(To)]^leq(^2)^2. = / #l^lH ^leqO^l)^!, (41) 



dr 
with which we can express g(r) formally as a Taylor series 



«T) = t^- <«) 



n n\ dr r 

All quantities (41) can be calculated analytically: 

'3\ J„(n\ Ai. 



m. „ 0.6760. m. = _ x .*m 



m = <^> = ^ « 0.6760. -^ = -1, -^ = 3,(0), 

d 3 o(0) d 4 o(0) 

On the other hand, the series (42) is not useful for the calculation of g(r) for a finite r since 
it is only an asymptotic one and it has the zero radius of convergence. So, to find the full 
autocorrelation function we resort to a numerical calculation. To be specific, we numerically 
solve the Fokker-Planck equation (39) for a new function S, 



9^(<Pi,t) 



dr 
with the initial condition 



IV-(£i,r), ^(<Pi,t) = / ^ 2 n[^ 1 (r)|^ 2 (0)]^ leq (^ 2 )^ 2 , (44) 
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^(<fil,0) = Pleq(<Pl)<Pl- (45) 

Then we can find g(r) from the integral 

flf(r) = / g?(^i~(^i,t). (46) 



The result is shown in Fig. 2. We find that g(r) asymptotically varies as g(r) oc e - 1 - 36859T . 
As expected, the numerical factor in the denominator is related to the lowest nonvanishing 
eigenvalue of (15): 



Ai « 1.36859W-^-#. (47) 

V 2471"^ 

Of course, g(r) is not purely exponential as in the A = case considered above, eigenvalues 
A n with n > 2 make some contribution to its expansion of the form (37). However, this 
additional contribution is small already at r = (e.g., \g~ 1 (0)-^-\ is only 8% more than its 
asymptotical value 1.36859 for r ^> 1), and it dies quickly because higher eigenvalues are 

significantly larger (in particular, A 2 ~ 4. 4537a/ -^^ H). As a result, g(r) becomes purely 
exponential with the error less than 1% for r > 0.3. £(t) = 0.5 at r ftJ 0.495, so the 
correlation time is t c ~ 7.62/ Hy\. 

C. Interacting case with two potential wells 

Another interesting case takes place when m 2 < and the potential has the form of a 
double well: 



J(* 2 - 


- <^o) 2 5 


^0 = 


1 ? 1 

\m \ 


A 



V^) = -(r-ti), ^ = -^- ( 48 ) 

In this model, the scalar held fluctuates around <f> = <f> in some domains and around <f> = — <f> 
in others. The correlation time, which is characterized by the lowest nonzero eigenvalue, 
Ai, of eq. (15), is much larger than in the massless A</> 4 theory. Let us assume y\H 2 <^ 
\m 2 \ <^ H 2 to demonstrate it. The effective potential in the Schrodinger equation (15) for 
the double- well potential (48) (made dimensionless using the quantities (38)) has the form: 

W{$) = Q~ 2 W(v) = ^ V + 7 ) 2 " \(^ 2 + 7), 



m 2 Q 2 87r 2 m 2 

^= —=mH- 2<0 - (49) 

It is depicted in Fig. 1 for 7 = —4. 

If the wells are sufficiently far from each other, then each energy level existing in one well 
in the absence of the other (including the lowest one in which we are interested in) becomes 
splitted into two with an energy difference between them being exponentially small. So, the 
eigenfunctions corresponding to A = and Ai may be represented, respectively, as 
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$oM = ^[$*M + $*H>)], (50) 

*iM = ^[^)-M-^)]. (51) 

where $*(<£>) is the normalized ground-state function for the single well localized around 
(f = (f b sa <f) with (^6 minimizing the potential W(ip). The problem we want to solve is close 
to that investigated in [19] (§ 50, Problem 3) in the WKB approximation. However, if we 
pretend to obtain not only the correct exponent, but the coefficient of the exponential, too, 
we may not directly apply the formula given there because the WKB approximation does not 
work for the lowest energy level generally. So, instead of the WKB approximation, we shall 
use the perturbation theory with respect to the exponentially small quantity S\ = 47r 2 Ai / H 3 . 
Two independent solutions of Eq. (15) in the static case A n = are: 

$(0) = e -«( V ) j $ (1) = e -v( V ) f V e 2v( V1 ) d(f ^ ( 52 ) 



The former function is actually proportional to that in Eq. (19). The corresponding static 
Green function of Eq. (15) is: 

G«(<p, Vl ) = sffMSfftei) - *S\<Pi)*$\<p). (53) 

Thus, the unnormalized eigenfunction for hrst non-zero energy level has the following form 
in the region ip ^> i7 2 |m| _1 (where the correction to the unperturbed wave function is small): 



$x {ip) = e~ v M - 2 £l / G st (<p, Vl )e~ v ^ d Vl 

J CO 

= e-< v Ul-2e 1 e 2v ^d Vl / e~ 2v ^ d^ 2 + 2e x / e ~ 2v ^ d Vl / e 2v( ^ d^ 2 . (54) 

Here the lower integration limit in the formula (52) for <& st is chosen to be <p = </> , i.e., to 
lie in the minimum of the potential V{cf>). However, any other value satisfying the condition 
(p 3> i7 2 |m| _1 may be used instead of it within the accuracy required to obtain a correct 
asymptotic value of A x . 

For (p <C </>o, the expression (54) takes the form: 

$i(v?) « e~ v( ^ (l + £l N l' V e 2v( ^ dtp^j . (55) 

The normalization coefficient N is defined in Eq. (19). In our case it is equal to: 



oo 



N= /"e-^)^ = J^n (56) 

V o^ \m\ 



oo 



(the main contribution is made by vicinities of the points ±</> )- Now $i(<£>) in Eq. (55) is 
represented as a linear combination of the two exact unperturbed solutions (52). Therefore, 
we may use this expression even if the second term inside the round brackets in its right-hand 
side becomes comparable to the hrst term. 
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On the other hand, $i(^) should be an odd function of tp } see (51). Therefore, $i(0) = 0. 
This gives the required equation for S\\ 

£l N [ °e 2t(B) ^ = 1. (57) 

Jo 

We, therefore, find 

A '= «T exp H^)' (58) 

The correlation time is t c ss In 2 A^ 1 that is exponentially larger than in the case of the 
massless A</> 4 theory. 

Similar to the massive non-interacting case, higher eigenvalues are given by the expression 

A 2ra , 2n+i ~ ^ rc n = l,2... (59) 

neglecting exponentially small splitting (the effective mass at the bottom of each well is 
equal to \/2 \m\). Using the spectral representation (37) and the formula (33), we arrive at 
the following expression for the autocorrelation function: 

i<- i 3# 4 / 2|m 2 |A , , 

TT 

.. exp(-A 1 t), t->-— (61) 

A \m z \ 

Eq. (60) shows the change in G(t) due to decay of correlations inside one domain, while Eq. 
(61) describes how the correlations disappear completely at very large times due to quantum 
jumps between different domains. 

Note that the exponent in Eq. (58) is nothing but the action for the Hawking-Moss (or 
de Sitter) instanton [20]. As discussed in [12], the physical sense of this instanton in the 
case under consideration becomes more transparent if one writes it in the static, "thermal" 
form: 

J S 2 = (1 _ ^2 r2) ^ + (1 _ H 2 r 2yl ^2 + f .2 ^2^ 

<f> = 0, (62) 

where t st is periodic with the period 2ttH~ 1 . Then it describes the "thermal" nucleation 
of a bubble of the phase with <f> = — <f> inside the phase <f> = <f> (or vice versa) through a 
quantum jump over the potential barrier which maximum is located at <f> = 0. The size of 
the bubble is equal to the size of the de Sitter event horizon H~ x . 

However, deriving Eqs. (58,61) we have got more — we have obtained the correct coeffi- 
cient of the exponential, too, that is equivalent to the calculation of a one-loop correction to 
the instanton contribution to G(t). This shows that not only does the stochastic approach 
go beyond the Hartree-Fock approximation, but it is able to reproduce results obtained using 
the instanton approach in the path integral formalism. 
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D. General case 

Here we present formulas valid for the case of an arbitrary potential V{cf>). Using the 
expressions (24-29) and repeating the calculation for the massless case presented above, 
we arrive at the following representation for the autocorrelation function in the equilibrium 
state: 



G(t) = j<pZ(<p,t)d<p, (63) 

where the function 

~(<p,t) = |^ 1/0leg (^)n [^(t)|^(0)] d 9l (64) 

satisfies the Fokker-Planck equation (12) with the initial condition 

Z(<p,0) = <pp leq (<p). (65) 

It is straightforward to generalize these formulas to the autocorrelation function of an 
arbitrary function F(ip). The result is that in order to obtain the corresponding expressions 
for the correlation function 

Gpdh-h]) = (FtffrhyFtffrh)]), (66) 

one has to substitute the multiplier ip in the right-hand sides of Eqs. (63-65) by F(ip). 
In other words, one has to solve the same Fokker-Planck equation with a different initial 
condition. 

Also useful is the spectral representation of G(t). For the initial condition (65), the 
coefficients (17) of expansion of £(<£>,£) into the complete set of eigenfunctions of Eq. (15) 
(which we denote by A n for this concrete case) have the form: 

A n = f ipp leq (ipy^<S> n (ip) dip = N- 1 f ipe- v ^<S> n (ip) dip, (67) 



where the normalization coefficient N is given in (19) and v is defined in (16). Note that 
A = if V(—ip) = V(ip). Then, using Eqs. (14) and (63), we arrive at the formula: 

G(t) = £ A n e- A ^ J <pe- v M$ n (<p) d^ = Nj2 ^e _A »' . (68) 

The relations (43) for the massless theory may be generalized, too. After differentiation 
of Eq. (63) and application of the Fokker-Planck equation, successive time derivatives of 
G(t) at t = follow: 



dt 2 24tt 2 AJ dip 2 T 24tt 

and so on (in the last equality in Eq. (70), we have inserted V(ip) = -| — \- m 2 y ). Another 
quantity for which an analytical expression may be obtained is I = J^° G(t) dt. Applying 
the method described in [11,12], we get 
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^fSf^MZ"^- 2 *' 1 ^) 2 (7l » 

(the expression is written for the case of a symmetric potential V(—ip) = V(ip) for simplicity). 
Eq. (71) may be used, e.g., to derive the formula (58) in a completely different way. 

A note should be made about a region of applicability of the expression for G(t). Ac- 
cording to the derivation of the Fokker-Planck equation in Sec. 3, we consider time intervals 
larger than i7 _1 |lne| ^> H~ x . So, strictly speaking, all formulas for G(t) are valid under 
this condition, too. However, G(t) is practically constant and equal to its value at zero lag 
G(0) = ((f) 2 ) for Ht <C min(i7 2 |m~ 2 |, A -1 ' 2 ) already. The expressions (69,70) for the time 
derivatives at zero lag should be understood in the same sense, i.e., actually they refer to 
the argument lying in the interval H -1 <C t <C H -1 min(i7 2 |m~ 2 |, A -1 ' 2 ). 

If A remains fixed and m 2 decreases from positive to negative values, two main trends 
follow from the results presented above: 

1) the autocorrelation function (in particular, G{0)) grows, 

2) the correlation time t c grows very quickly. 

To visualize these trends we made a numerical calculation of the dimensionless correlation 
function g(r) = Q 2 G(t) for the cases 7 = 1,0,-1,-4 using the dimensionless variables 
defined in (38,49). The results are depicted in Figs. 2 and 3. They confirm the expected 
behaviour. The case 7 = —4 is already close to the asymptotic case considered in subsection 
C, in particular, the lowest non-zero eigenvalue Ai differs by less than 10% from that given 
by Eq. (58). 

V. GENERAL EQUILIBRIUM TWO-POINT CORRELATION FUNCTION AND 

ITS DE SITTER INVARIANCE 

Let us now consider the spatial correlation function 

G(r,t) = (t(x 1 ,t)t(x 2 ,t)), r = | Xl -x 2 |. (72) 

Repeating the same procedure to derive (12), one can also find the Fokker-Planck equation 
for the two-point PDF at equal time / o 2 [ ( ^'i(xi, t), (,c 2 (x 2 , t)] as 

f)n /7 3 r) 2 n 

-7^i(xi,t),^ 2 (x 2 ,t)] = T vl p 2 + T V2 p 2 +j (ea(t)H\x 1 -x 2 |) — 2 . (73) 

Since jo(z) is a rapidly oscillating function for \z\> 1 and we are concerned with coarse- 
grained quantities, we can, and in fact we should, adopt the following approximation: 

joMt)#| Xl - x 2 |) ~ 0(1 - ea{t)H\^ - x 2 |), (74) 

which is employed hereafter. The above replacement corresponds to the approximation that 
two points Xi and x 2 with their proper distance a(t)|xi — x 2 | = a(t)r < (ei7) _1 have 100% 
correlation while those with larger separation are mutually independent at the time t. 
In fact, one can readily see that 

P2(<Pl,<P2) = Plq(<Pl)6( l Pl- l P2) (75) 
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constitutes a static solution of (73) with j replaced by unity. Therefore one can find the 
equilibrium two-point PDF with a(t)r > (ei7) _1 by solving 



dp 2 
dt 



[(^ 1 (x 1 ,t),(^ 2 (x 2 ,t)] = T vl p 2 + T V2 p 2 , (76) 



with the initial condition 

P2[<f 1(^1, U),<f 2(^2, U)] = Piq(<Pi)t>(<Pi ~ ^2), Ht r = -ln(ra eH). (77) 

Now we turn to the calculation of the general two-point correlation function 

G(r } t 1} t 2 ) = (^(x!,*!)^,^)). (78) 

Following essentially the same line of arguments, we arrive at the basic expression for the 
general two-point PDF in the stochastic approach: 



/ 2 [^l(Xl,ti),^2(x2,t 2 )] = / n[^i(Xi,ti)|^ r (Xi,t r )]n[^2(x2,t 2 )|^ r (x2,t r )] / Oi q (^ r )^ r , 

t r = — H~ x ln(ra eH), r = |xi — x 2 |, (79) 

where II[(^i(x, ti)\tp 2 (x. } t 2 )] satishes the Fokker-Planck equations (26,27) with respect to 
both its time arguments and the initial condition (28). Note also that the spatial points Xi 
and x 2 were inside the same elementary averaging volume at the moment t = t r , thus, they 
had the same one-point PDF at this time. 

Of course, it is assumed here that t 1? t 2 > t r . This means that the formula (79) refers to 
space-time points (xi,ti) and (x 2 ,t 2 ) lying outside each other's light cones. On the contrary, 
if these points may be connected by a time- like or null geodesies (z > f in Eq. (35)), then 
min(ti,t 2 ) < t r that corresponds to the points (xi, min(ti, t 2 )) and (x 2 , min(ti, t 2 )) lying 
inside one elementary averaging volume, i.e., coinciding in terms of the stochastic approach. 
Then the 2-point function G(r } ti 7 t 2 ) is equal to the autocorrelation function G(\ti — t 2 \) 
found in the previous section. 

Returning to the case of a space-like separation between two points and using the spectral 
decomposition (14) for Tl[(pi(x. } ti)\(p 2 (x. } t 2 )] together with the definition (67) of A n , we 
obtain 



G(r } t 1} t 2 ) = / diftift / d<p2<P2P2[<Pi(x-i,t 1 ),<p 2 (x2,t 2 )] 

= jj d^ d^ 2 ^^er^er^) J d^ r e 2v ^p lci (^ r ) 

X £ $ m (^)<M^)e- Am(fl - M E M^)^ n ^r)e- A ^- U) 

ra n 

= ^EE /%$ m (^)$n(^)^A.e" Ara(tl " ir) e" A ' ,(f2 " (r) 

m n 

= NJ2 A 2 n e- A ^ tl+t2 - 2t ^ = NJ2 A 2 n e- A ^ tl+t ^ exp (-2A n #- x ln(ra eH)) . (80) 



Now we may choose the small parameter e of the coarse-graining (6) in such a way that 
e -2A n /H r^ ^ £ or a rj A n which make a significant contribution to G(t). This is achieved if e 
satishes the inequalities 
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exp (-min(# 2 |m- 2 |,A- 1/2 )) < e < 1. (81) 

As a result, we may omit e Irom the last line in Eq. (80). Then, from the comparison of (80) 
with (68), the expression for the general two-point correlation function in the equilibrium 
state in terms of the autocorrelation function G(t) found in the previous section follows: 

G(r, t u t 2 ) = G (h + t 2 + 2H- 1 ln{ra H)) . (82) 



Using the de Sitter invariant form z(xi 7 x 2 ) (35) in the regime of a strongly space-like sepa- 
ration between two points (z < 0, \z\ ^> 1) to which the formula (79) refers, Eq. (82) can 
be rewritten in the form 

G(r,t 1 ,t 2 ) = G(H- 1 ]n\2z-l\) (83) 



that applies both to the cases of large space-like and time-like separations. As explained at 
the end of the previous section, we may use this expression even in the vicinity of the point 
z = 1 where the right-hand side of (83) is simply equal to G{0) = (<f> 2 )- Therefore, we have 
proved the de Sitter invariance of the two-point correlation function. 
In particular, the equal-time spatial correlation function (72) 

G(r,t) = G(r,t,t) = NJ2 A 2 e- 2A " f exp (-2A n i7" 1 ln{ra H)) = A]T A 2 n (RH)- 2A »' H (84) 

n n 

appears to depend on the physical spatial distance R = ra e Ht only. Thus, a typical scale of 
spatial correlations in the equilibrium state, which defines the characteristic size of domains 
with positive and negative ip and domain walls with ip ~ between them, does not expand 
with the expansion of the t = const hypersurface of the de Sitter space-time. Similar to 
the correlation time t C} we can introduce the spatial physical correlation radius R c using the 
definition G(R C ) = G(0)/2. Then the de Sitter invariance property (82) gives the relation 
between t c and R c : 

R c = H- 1 exp ( ^) . (85) 



VI. DISCUSSION 

In the present paper, we have further developed the stochastic approach to inflation by 
constructing a method of calculation of the 2-point and higher-point correlation functions 
of a quantum scalar held in the de Sitter background in terms of solutions of the Fokker- 
Planck equation (12). This enables us to demonstrate the existence of the equilibrium state, 
which has the de Sitter invariance, for a self-interacting scalar held with a small mass term 
of an arbitrary sign. It is easy to see that this state is the unique attractor for any initial 
configuration. For example, the temporal two-point PDF p 2 [tpi(x. } t 2 + t) } tp 2 (x. } t 2 )] with 
an arbitrary initial distribution / 0i[( J c 2 (x, t 2 )] is obtained solving the Fokker-Planck equation 
(29) with the initial condition 

/ o 2 [^i(x,t 2 ),(^ 2 (x,t 2 )] = / o 1 [(^ 2 (x,t 2 )]<5((^ 1 - ip 2 ). (86) 
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Since any solution p 1 [cp 2 {x,t 2 )] approaches piq^) as t 2 — > oo, / o 2 [ ( ^i(x,t 2 + t),cp 2 (x,t 2 )], 
as well as the autocorrelation function G(t), approaches the equilibrium counterpart in this 
limit. Moreover, as long as we are interested in a fixed physical length scale, similar argu- 
ments hold for the spatial correlation function as well, and the asymptotic spatial correlation 
function does not depend on time even though the space (i.e., the hypersurface t = const in 
the system of reference (3)) is expanding. It is evident that had we sticked to the conven- 
tional perturbative expansion around the free asymptotic state, we could not have obtained 
the above results. 

Thus, any weakly interacting scalar held with a small mass eventually loses the memory 
of its initial state at the beginning of the de Sitter era as t — > oo. The time scale for this to 
happen may be characterized by a relaxation time t re \. Its value depends on the quantity 
involved. Generally, the relaxation time is of the order or less than the correlation time t C} 
e.g., t re i ~ t c ~ In2/Ai for the two-point PDF. On the other hand, repeating the calculations 
in Eqs. (68, 80) for an arbitrary non-equilibrium initial one-point PDF pi(tp } t ) } we arrive 
at the result that the relaxation time for the two-point correlation function t Te \ ~ A^ 1 for 
any reflection symmetric potential V{—cf>) = V{cf>). So, in this case 

t rel ~mm(H\m- 2 \ } H- 1 \- 1/2 ') , (87) 

and t re i may be much less than t C} cf. Eqs. (58) and (59). 

Very important is the result (84) showing that the equal-time spatial correlation function 
depends on the physical radius R only (and not on the comoving radius r). This has dramatic 
consequences for the general picture of fluctuations of the scalar field which are especially 
interesting in the case of a potential with two wells and \m 2 \ ^> y\H 2 (Sec. 4, subsection 
C). In the latter case, the space is covered by clearly pronounced domains with ip ss ±</> 
with relatively thinner (but much thicker than H~ x ) domain walls between them. Typical 
size of the domains in the equilibrium state is given by the spatial correlation radius (85). 
It is fantastically large because it has the form of a double exponent, see Eq. (58), but 
still it is finite. The size of the domain walls Rd w is determined by A^ and contains single 
large exponent: \n(Rd w H) ~ H 2 \m~ 2 \. Both these physical sizes do not depend on time on 
average that corresponds to their typical comoving radius r shrinking with time. Therefore, 
the expansion of the background physical space volume oc e 3Ht is not accompanied by the 
corresponding expansion of individual domains, instead of it, new domains separated by new 
domain walls are created with constant rate. This purely quantum (or stochastic) behaviour 
has to be contrasted with the deterministic classical expansion of inflationary domain walls 
well studied in the context of the new inflationary scenario (recently, renewed interest in 
this problem arose, especially, in connection with more complicated topological defects like 
strings and monopoles [21,22]). Though our calculations refer to the case of a test field 
in the stable de Sitter background, it is clear that one may expect a similar effect when a 
quasi-de Sitter background is produced by the scalar field itself, namely, that the classical 
monotonous quasi-exponential expansion of an individual domain wall at the first stage of 
inflation stops at the late asymptotic stage due to the loss of coherence produced by the 
quantum noise, and new domain walls are created instead of it. This will be considered in 
more details in another publication. 

In this paper, we assumed the exact, "eternal" de Sitter background with H = const. In 
real inflationary models, H is not exactly constant. The characteristic time of its change is 
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Tjj = H/\H\ ^> H~ x . Then the behaviour of a quantum self- interacting scalar held in such 
a background depends on the relation between tjj and the relaxation time. If t re \ <C tjj, the 
scalar held has enough time to reach the equilibrium state described in this paper for each 
successive Hit). In the opposite case, evolution of the scalar held is significantly different 
and constitutes a separate problem. 
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Figure Captions 

Figure 1 Dimensionless effective potential W[tp] for a quartic potential V[(f)] = -^(f) 4 (solid 
line), and that for a double well potential V[<f>] = |(</> 2 — <f>l) 2 = ^(<f> 2 + ~fQ~ 2 ) 2 with 
the dimentionless mass parameter 7 = —4 (dashed line). 

Figure 2 Dimensionless autocorrelation function g(r) of a scalar held with a potential 
■^(4> 2 + ^Q~ 2 ) 2 . Three curves represent the cases with 7 = —1,0, and 1, respectively, 
from the above. 

Figure 3 Same as Figure 2 but with 7 = —4. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9407016vl 



25 

W«>) 
20 



15 



10 







-5 



A 

_ I 

I 

~~ I 

— I 

- I 
I 

I I 
r i 
■j- i 

I I 
T i 












i 

i 

i 

i 

i 

i 

i 

i 

i 

i 


— 




i 
4 i 

I i 
J i 

| i 
— i 

— L i 

_ i 

i 

I i 

— \ 

\ 










i 

i 

i 

i 

i 

i 

i 
i 
i 
i 


— 












i 


I — 


\ 










i 






\ 








i 






\ 








i 




— \ 


\ 
\ 








/ 


/ — 


— \ 


\ 








i 


/ — 


— \ 


\ /^ 


^^^^ 






i 


/ — 
















\ / 






\/ 
















— 


— . 










/ \. ^- 


— " 










— \ 












— 


: 














: 



-2 




Figure 1 



* 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9407016vl 



1.2 



1 



g(x) 



.8 



.6 



.4 



.2 












1 



2 
Figure 2 



This figure "figl-3.png" is available in "png" format from: 



http://arXiv.org/ps/astro-ph/9407016vl 



4.2 



g(x) 



3.8 



3.6 



3.4 



3.2 








1 



2 
Figure 3 



